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Enrichment of nuclear spin isomers of molecules by infrared radiation resonant to 
molecular rovibrational transition is considered. Special attention is given to the 
enrichment by light-induced crossing of far separated ortho and para states. 

INTRODUCTION 

Nuclear spin isomers were discovered in the late 1920s when ortho and para H2 were 
separated by deep cooling of hydrogen gas [|l| . This technique was successful because of 
extremely large gap between the ortho and para states in H2. Investigations of heavier 
isomers were not possible to perform until recently because of the lack of appropriate 
enrichment methods. First separation of heavy molecules (CH3F) was achieved using 
the Light-Induce Drift effect §. Later, the isomers of H2O §, H2CO [§, U2 §, and 
C2H4 were enriched. A few practical applications of spin isomer are already known. 
For example, isomers can affect chemical reactions |]^, can significantly enhance NMR 
signals |T0[, or can be used as spin labels in NMR based techniques. 

Further investigations of nuclear spin isomers and progress in their applications need 
new enrichment methods. This is rather difficult and challenging task because physical 
and chemical properties of spin isomers are almost identical. Recently new approach 
to the isomer enrichment based on molecular excitation by strong electromagnetic field 



was proposed [11, 12, [T^, 14 1. In this paper we discuss the isomer enrichment based on 
the coherent control by strong infrared radiation. We make emphasis on the conditions 
necessary for the experimental observation of the phenomenon. 
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LEVEL SCHEME 

The level scheme is shown in Fig. |l]. We assume that test molecules have only two 
nuclear spin states, ortho and para. This is very often the case, e.g., in CH3F ITSl, p!6|. Each 
of the two vibrational states (ground and excited) has one ortho-para level pair mixed 
by intramolecular perturbation. Molecules are embedded into a "nonmagnetic" gaseous 
environment characterized by zero cross-section for collisional ortho-para transition. The 
infrared radiation is resonant to the rovibrational transition m — noi the ortho molecules. 

The molecular Hamiltonian in our problem consists of four terms, 

H = Ho + hHrad + fiV + hV'. (1) 

The main part, Hq, has the eigen ortho and para states shown in Fig. |1[ hHrad describes 
the molecular interactions with the external radiation that will be taken in the form of 
monochromatic travelling wave, 

Hrad = -2G cos{uLt - kr) ; G = Eod/2;i, (2) 

where Eq; ujl and k are the amplitude, frequency and wave vector of the electromagnetic 
wave, respectively; d is the molecular electric dipole moment operator. V and V' are the 
intramolecular perturbations that mix the ortho and para states in excited and ground 
vibrational states, respectively. These perturbations originate from the nuclear spin-spin 
and spin-rotation interactions Jl^, |18|, |19], pOl . 

Suppose that a test molecule is placed initially in the ortho subspace of the ground 
vibrational state. Due to collisions the molecule will undergo fast rotational and velocity 
relaxation inside the ortho subspace. This will proceed until the molecule jumps to the 
state m', which is mixed with the para state k' by the intramolecular perturbation V' (see 
Fig. p, or to the state n which is mixed with the para state k by the combined action of 
the external field and intramolecular perturbation V. Admixture of a para state implies 
that the next collision can move the molecule to another para states and thus localize it 
inside the para subspace. The spin conversion occurs. 

The molecular center-of-mass motion affects the process because of the Doppler effect. 
Vibrational degree of freedom plays an important role too. After ortho-para transition in 
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upper vibrational state the molecule relaxes to the ground vibrational state. Consequently, 
back conversion can only occur through the unperturbed m' — k' channel. It is essential 
that conversion through the m—k channel can be substantially affected by electromagnetic 
radiation. 

Quantitative description of the isomer conversion can be done using the perturbation 
theory The zeroth-order theory neglects the intramolecular ortho-para state mixing 
by V and V' and thus contains no ortho-para conversion. In this approximation an 
external field does not affect the population of para states. The populations of ortho 
states in the model of strong collisions read, 

Po(e,J,v) = Po{p{fl)TevWjf{\) +p{Q,Y)[TetWj + TeSjm]}, 

PoigJ,^) = PoWjfiy) - Po{p{VL)Tg^Wjf{y) +p{VL,^f)[TgtWj + Tg5jn]}- (3) 

Here po is the total concentration of ortho molecules; e, g designate the excited and ground 
vibrational states, respectively; j is the rotational quantum number; Wj is the Boltzmann 
distribution; /(v) is the Maxwell distribution and v is the molecular center-of-mass ve- 
locity. The level populations in Eq. are split into the terms nonequilibrium in j and 
V, nonequilibrium only in v, and equilibrium in j and v. The magnitudes of these terms 
are proportional to specific relaxation parameters r's (see Ref. for more details). The 
function v) is the excitation probability of molecules with particular velocity v and 
p(f^) = /p(fi,v)rfv. 



Y{n) 



r| + (i]-kv)2' 



Ti 1 + K 1 + Y{n)T2T^ r (iZ) + [il - kv)^ 

Here Q = —uj and u = Umk is the gap between the states m and k. The function p{Q) is 
the Foigt profile of the absorption line having the homogeneous linewidth Tb = + k, 
where k is the saturation parameter, k = 2|G'prir~^; \G\ = i?o|'imn|/2/i is the Rabi 
frequency; T is the collisional linewidth of the absorbing transition; ri = TetWm + TgtWn + 

Te + Tg and T2 = TeyWm + TgyWn- 

In addition to the level population change, an external radiation creates a coherence 
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(off-diagonal density matrix element) on the m — n transition, 

p,(n|m; v) = -^Po^^[r - - kv)]e'(™-'^'-). (5) 



ENRICHMENT 

The first-order theory takes the intramolecular perturbations V and V' into account. 
In this approximation the spin conversion appears and it can be affected by the radiation. 
An external radiation changes the concentration of ortho molecules from the equilibrium 
value, Po(0), to the steady-state value, p^. Because the total concentration of the test 



molecules is conserved, one has steady-state enrichment of para molecules |[21 



/? = - 1 = 1 - 7 = 7/ree - 7n + 7n + 7coh. (6) 

Pp(0) 7 

Here an equal concentrations of ortho and para molecules in equilibrium gas has been 
assumed that, implies Wm' = wy- The conversion rate 7 has four terms of different origin. 
The first one, 

2V\V' 



7 free = , ,^ [Wk' + Wm'], (7) 



is the conversion rate without an external electromagnetic field; u' = uJm'k' is the gap 
between the states m' and k'. One can see from Eq. that the light-induced enrichment 
is enhanced by small 7/ree- The second term is due to the light-induced level depopulation 
in the ground vibrational state, 

2r|i^'p 1 

The last two terms in 7 is convenient to present in a form that can be easily obtained 



from pTI. The "noncoherent" contribution to the rate. 



L'^+UJ^ V L + too / J 1 1 + z(iZ + CJi — kv) 

Here r = TevWm + TetWm + Te- The "coherent" part, ■jcoh, originates from the light-induced 
coherence (^, 

\V\M^) , \VW T + T,+tu, f p{n,v)dv 
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The new homogeneous width, Fi, and the ortho-para level gap, ui, are 

One can see from Eqs. (^), (|^), and (p!0| ) that the conversion rate, 7, has two peaks. The 
peak at the absorption line center (f2 = 0) appears because the excitation probabilities, 
p{fl) and p{fl, v), have maximum at f2 = 0. The homogeneous width of this peak is equal 
to that can be very large in molecular systems due to the slow rotational, velocity, 
and vibrational relaxations. The peak at Q = —uoi is due to the light-induced crossing 
of m and k states (see Ref. [l^, ^) and appears formally because of the resonant 
denominators in the integrands in Eq. (P) and (|T0|) . The homogeneous width of the peak 
at f2 = —uji is equal to Fi. 

Coherent control allows to achieve substantial enrichment if the mixing in ground 



and vibrationally excited states has the same order of magnitude On the other 

hand, close coincidences between ortho and para states in molecules are rather rare. 
Consequently, such close level pairs in excited vibrational state may be absent, or they 
cannot be accessible by the available radiation. Below we estimate the enrichment in the 
case of large ortho-para gap in excited state, uj ^ kvQ ^ F, where kvQ is the Doppler 



parameter and vq is the molecular thermal velocity, vq = y2kBT/M. For moderate 
radiation intensities, \G\ <^ cu, the field-dependent terms at the absorption line center, 
r2 = 0, are 

7;(0) ~ 0.57/.eeP(0)(v + rgt); 7n(0) ^ ^p(0)rF; 7,„,(0) ^ -^P(O). (12) 

The terms 7n(0) and 7co/i(0) have the "off-resonant" nature which manifests itself by the 
large ortho-para gap, u, in the denominators. At large u, the conversion rate, 7(0), is 
determined mainly by 7^(0) and the radiation depletes the para molecules, 

/3(0)~-0.5p(0)(v + r,i). (13) 

This is analogous to the enrichment by light-induced population change |[Tl|, |12[. Maxi- 
mum of /3(0) is achieved at sufficiently strong field, F^ ^ kvo- In this case, Y{0) ^ 1, 
p[0) ~ Wn/ (ti + T2) and enrichment reaches the value, 

(3^U0) ^ -wj-^i^^^. (14) 

Tl + T2 
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Depending on the relaxation scheme and the origin of the relaxation "bottle neck", one 
can achieve different level of enrichment. In the most favorable situation of fast velocity 
and rotational relaxations the enrichment reaches -1/4. 

Larger enrichment can be achieved at the level-crossing resonance. The field-dependent 
rates at the Q = —Ui can be readily estimated if Tb ^ kvQ. In this case p(— Ci;i,v) ~ 
p(— co'i)/(v) and one has, 

7„(-c^i) ~ ^, ' ' ' , p{-ooi)t, 7coh(-^i) ^ ^ ' ' p{-u,)r-\ (15) 

The terms, 7„(— c^i) and 'jcohi—^i) are due to the light-induced crossing of the m and 
k states that has an efficient width equal kvo/\^ (compare with the level crossing in 



molecules being at rest [Q). The main contribution to 7(— c^i) comes from '~fcoh{—^i) if 

The results (|12D and (|15D allows to determine which peak, Q = 0, ot Q = —u, should 
be used for enrichment. At low radiation intensity larger enrichment is achieved at Q = 
because p(0) ^ p{—uJi) if ^ u;. If radiation intensity is high, ~ cj, larger 
enrichment is obtained a.t Q = —uji. The ratio of the two peaks reads, 

2(z^^24l^.V^r-V,. + r„)-. (16) 

7(0) 7n(0) IfreekVo 

The amplitude of the enrichment peak at Q = —uji as a function of IGj and various 
values of uj are shown in Fig. The calculations were done using the exact expressions 
for the conversion rates, (^-([l^)- In this numerical example the following parameters 



were used V = 3 kHz, V = 5 kHz, F = 2 MHz, oo' = 130 MHz, kvo = 30 MHz; 
Ti = 5 MHz~^; T2 = 1 MHz^^ and all Boltzmann factors equal, Wm = Wn = Wk = Wm' = 
Wk' = 10^2_ 



DISCUSSION 



Conditions necessary for the experimental observation of the light-induced enrichment 
are not easy to fulfil. First of all, the field free conversion rate, 7/ree; of the molecules 
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under consideration should be low in order to enhance the role of light-induced changes 
in the total rate, 7. Note that the rate ^free should combine all significant contributions, 
e.g., caused by the molecular collisions with the walls of the gas container. At present, 
only in a few rare cases (see the Introduction) the rates 7/ree are known. 

Next task is to determine the target ortho-para level pair in the excited vibrational 
state. For this, one has to know rather accurately the relative positions of ortho and para 
states and the intramolecular ortho-para mixing efficiency as well. 

Perhaps, the most difficult problem is to find powerful radiation having right frequency. 
The most promising sources of such radiation are the powerful Free Electron Lasers. 
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FIG. 1: Level scheme. Vertical solid line indicates infrared excitation. Vibrational relaxation is 
shown by dashed vertical lines. V and V indicate the intramolecular ortho-para state mixing. 
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FIG. 2: Enrichment of para molecules as a function of Rabi frequency, |G|, and various ortho- 
para level gaps in excited vibrational state, lo. 1) lo = 100 MHz, 2) 500 MHz, 3) 1000 MHz, 4) 
1500 MHz, 5) 2000 MHz. 



